In-medium modification of P-wave charmonia from QCD sum rules 
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We investigate the changes of the masses and widths of XcO and xa in hot gluonic matter near 
Tc and in nuclear medium using QCD sum rules. As in the previous works for the J/ip and r^c, 
in-medium effects are incorporated through the changes of gluon condensates. Twist-2 terms for 
the 'Po(xco) and ^Pi{x ci) are also included for the first time. The results show that larger mass 
shifts and width broadenings take place as compared to the S-wave states. As the critical change 
take place near Tc, related measurements can reveal critical phenomenon in QCD. 



0^ 
O 
O 

a 



^ ■ 

Oh' 
I ■ 

Oh' 

CD 



(N 
> 

m 
in 



00 

o 

00 

O 



X 



PACS numbers: 24.85.-|-p,14.40.Gx,11.55.Hx,12.38.Mh 



I. INTRODUCTION 

Spectral property of charmonia in medium can provide 
experimental information about the confinement nature 
of QCD. At high temperature, hadronic matter under- 
goes a phase transition into the deconfined phase. Re- 
cent experimental data measured at Relativistic Heavy 
Ion Collider (RHIC) reveals that the deconfined matter 
is strongly interacting [E H, H, 01 • It was previously ex- 
pected that the production of J/ tjj will be suppressed by 
the color screening in the deconfined matter formed in 
the initial stages of heavy ion collisions [1, On the 
other hand, recent experimental observation shows non- 
trivial suppression pattern that are not easily understood 
by a simple suppression mechanism. An important fact 
in J/ip production is that they are not only produced 
directly but also through decay from higher mass states 
such as tp' and Xc- Since such excited states are likely to 
dissolve at just above the critical temperature Tc whereas 
the ground state does at 1.6Tc, as have been shown by lat- 
tice calculations 0, B , sequential suppression scenario 
[lo| has been proposed, in which the initial suppression 
comes from the disappearance of the feed down from ip' 
and Xc followed by the suppression of directly produced 
J/ip. In previous works, two of us (K.M and S.H.L) have 
studied the properties of J/ip and rjc using QCD sum 
rules, and have showed that their masses and widths 
change following the critical behavior of gluon conden- 
sates, which are extracted from lattice QCD [Til. [T^. In 
this paper, we apply the same formalism to the Xco and 
Xci states to investigate how the properties of these states 
are affected by the QCD phase transition, and explore 
the consequences in the experimental measurements. We 
also study the nuclear matter case which may serve as a 
testing ground for the phase transition through precursor 
phenomena. 

The paper organized as follows; in the next section, we 
briefly summarize our formalism and present new Wilson 
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coefficients. In section lllll we present quantitative results 
of in- medium change of masses and widths of XcO and Xci ■ 
Section [IVI is devoted to discussion and summary. 



II. FORMALISM 



QCD sum rules are based on the current correlation 
function and the dispersion relation [l3|. In this paper, 
we consider the scalar current ~ cc and axial vec- 
tor current = {qt^qu/q^ - 9t,v)cYlbC for XcO and xa, 
respectively. For finite temperature we start with the 
correlation function 



Ii-'{q)=i d'xe^'i-^T]{x)m)l 



(1) 



with (• • • )t being the Gibbs average {0{x))t = 
TT{e~^^0)/TY{e~^^). For the scalar, there is only one 
independent polarization function, but for the axial vec- 
tor correlation, there are two independent polarization 
with respect to the medium. This is so because when the 
axial vector meson is moving with a finite three momen- 
tum q with respect to the medium at rest, the response 
will be different depending on whether the polarization is 
parallel or perpendicular to the three momentum. Here- 
after, we assume that both the medium and the cc pair 
are at rest so that q = {ui, 0). Then, for both the scalar 
and axial vector meson, there will be only one indepen- 
dent mode. For the axial vector correlation function, 
one can extract the scalar correlation function through 
^■^(g) ^ H^(g)/g2 and n^(g) = -Yl^" / {3q'^) . Taking 
q^ < 0, the dispersion relation relates the correlation 



function with the spectral density p(s) [U. 
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In the QCD sum rules, we calculate the correlation func- 
tion by means of the operator production expansion 
(OPE) with condensates while we model the hadronic 
spectral density with a pole and continuum. In the 
present calculation, we neglect the continuum part in the 
hadronic spectral density since we can suppress the con- 
tribution by means of the moment sum rules. Recently 
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the role of the zero mode in the spectral density has been 
emphasized in several literature [Tsl [iGj since it gives 
a constant contribution to the imaginary time correla- 
tor from which the spectral density is reconstructed via 
maximum entropy method in lattice QCD. This comes 
from Landau damping effect, i.e., scattering of the cur- 
rent with on-shell quarks in medium and has been known 
as the scattering term 17] in the framework of QCD sum 
rules. This term is also related to the transport coef- 
ficients as they are the terms responsible for the long 
range order surviving at zero momentum and energy of 
the current. The scattering term exists whenever the 
thermal excitation, thermal charm quark in our present 
case, couples directly to the current, and contributes to 
the spectral density at zero energy in the kinematical 
limit of q — !■ with a strength proportional to the ther- 
mal charm quark density at finite temperature p^ . It 
can be shown that the same scattering term can be ob- 
tained in the OPE side by resuming the thermal charm 
quark contribution to the zeroth order charm quark oper- 
ators, that did not have direct charm quark contribution 
at zero temperature as they were converted to gluonic 
operators in the heavy quark expansion figj. Therefore, 
while the scattering term contribution is non-negligible in 
the spectral density, similar term in the OPE will cancel 
its contribution and thus will not be considered further in 
this work. In a sense, it is the advantage of the QCD sum 
rule method to artificially turn off the effects of thermal 
charm quarks altogether and concentrate on the thermal 
gluonic effects on the pole. For confined phase, there is 
no on-shell quarks and we do not have to consider the 
scattering term from the beginning since we have only 
glueballs other than the heavy quarkonia in this phase. 
Therefore, we will neglect the scattering terms in both 
the confined and deconfined phase. 

For heavy quarkonia, the moment sum rule gives a 
systematic procedure. The moment is defined as 



d 
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(3) 



Noting that the expansion parameter in the OPE for 
heavy quarkonia is AQgj-,/(4rn^ -I- Q^) and focusing on 
temperature as large as Aqcd so that the expansion pa- 
rameter is not much modified from the vacuum value, 
we can truncate the OPE at the leading order pertur- 
bative correction and the dimension four giuon conden- 
sates, with the dimension four twist-2 gluon condensates 
being a new addition as compared to the vacuum case. 
In this approach, all the temperature effects in the OPE 
are encoded in the temperature dependence of operators. 
Therefore the moment ([3]) for the OPE side becomes, for 



Mn(e) = A.(e)[l + an(0as(0 + &«(0'^fc + cn(0'^c], (4) 



where An, Q!„, 5„ and c„ are the Wilson coefficients for 
the bare loop, perturbative radiative correction, scalar 
gluon condensate and twist-2 gluon condensate, respec- 
tively. In Eq. ([4]), 0h and 4>c have temperature depen- 
dencies through the gluon condensates. They have been 
extracted from lattice calculation of pure SU(3) theory 
and shown in Ref . [l3] ■ Here we would like to stress that 
they are not order parameters in a strict sense but their 
behavior more or less reflect the phase transition, i.e., 
sudden change across the critical temperature Tc. 

As for the Wilson coefficients. An, an and 6„ are listed 
in Ref. [2^ . c„ for the scalar and axial currents are cal- 
culated for the first time in this work by making use of 
the background field technique [llj. Defining the scalar 
functions 11^ (g) = n^(g)/g2 and n^(g) = for 
the scalar channel and the axial vector channel correla- 
tion function, respectively, we find that the correction 
term for the twist-2 gluon operators in OPE becomes 



Q2 g4 



(5) 
(6) 



where y = Q^/m? and Jn{y) — dx[l + x{l — x)y\ ". where p ~ and F{a,b,c; z) denotes the hypergeo- 
Therefore c„ is given by metric function 2Fi{a,b,c; z). One can see that c„ can 



c, 



An{n + 1) F{n + l,\,n+^;p) 



be expressed as a sum of 6„ and an additional term, as 



Sj-^) ^ ^S(-^-) ^ -t'H"--'- V'--r 2'"-^ 2'^^ ^ (j^ in the case of the vector and pseudoscalar case [22], and 

3(1 + F{n, 2 1 + 2 ' ^) that the additional terms of the scalar and axial vector 
-I- 

~^ ^) -^"("-7 2 ' "I" 2 ' Z') Masses and widths are extracted in the same manner 



^r^^ ^.Ar,^ 4n(n + 1) + 1, i, n |; p) channel are identical 
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as in the previous works for J/ip and 7}^ ■ Taking 

ratio of the moments and equating the OPE side with 
the phenomenological side, we find 



Mn-l 



OPE 



(9) 



phcn. 



In the phenomenological side, we simply express the pole 
term with the relativistic Breit-Wigner form and neglect 
the continuum contribution. Namely, we put 



1 



p{s) = - 



TT (s ■ 



sr2' 



(10) 



This simplification is justified by taking the moment 
since it suppresses the high energy contribution in the 
phenomenological side. While the continuum slightly 
modifies the mass value, its contribution does not af- 
fect the in-medium change from the vacuum value, which 
is the interest of this paper. Also, the thermal factor 
tanh(y^/2T) is put to unity since we are considering 
temperatures much lower than the mass of charmonium. 



III. RESULTS 
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FIG. 1: (Color online). Masses in the F — > limit as a 
function of n and T/Tc- The upper panel is for XcO and the 
lower for Xci. 

Following Ref. [1^, we work with ^ = 2.5 to validate 
the OPE since the Wilson coefficients are larger than 



those of the S'-wave. The parameters are the same as in 
Ref. [H, with as{^ = 1) = 0.21 and = 1) = 1-24 

GeV. Since we have neglected the continuum, vacuum 
masses are a little larger than the experimental data. 
However, it does not affect the results of medium- induced 
changes discussed in what follows, as we will discuss later. 

Figure [T] shows the masses in the F — > limit cal- 
culated from = M„_i/M„|ope — 4to^^ at various 
temperatures. At large n, the OPE breaks down, while 
at small n, the continuum contribution becomes impor- 
tant. Therefore, the region where both approximations 
are good enough and the sum rule is reliable shows up in 
the plot as a minimum, the value at which corresponds to 
the physical mass 20, 2^]. As seen in the S'-wave cases, 
the masses decrease as temperature increases. The sum 
rule can be applied as long as the moment ratio has a 
minimum. We can see that the OPE breaks down at 
T > l.OSTc for both XcO and Xci- To see the behav- 
ior clearly, we also display each term in Eq. (|4]) in Figs. [2] 
and[31 One can see that the OPE breaks down because of 
large expansion coefficients which should be smaller than, 
typically, 0.3. At the marginal temperature T = 1.03Tc, 
the perturbative correction term is close to 0.5; this is so 
because the minimum of the moment ratio is shifted to 
the large n direction where the correction terms increase 
as can be seen in Fig. [3l Therefore, the result at this 
temperature is less reliable than that at lower tempera- 
tures. From Fig. [31 one also notes the opposite sign be- 
tween bn4'b (scalar gluon condensate) and c„0c (twist-2). 
This is a general feature of the condensates contribution 
coming from the fact that the scalar condensate remains 
positive at temperatures not so higher than Tc while the 
twist-2 operator is always negative. Their magnitudes 
almost coincide at T = 1.03Tc and the twist-2 term ex- 
ceeds the scalar condensate at T = 1.04rc- The fact that 
the individual OPE term is large, and that the pertur- 
bative as contribution becomes relatively larger than the 
scalar condensate, lead to the loss of stability of the OPE 
at the marginal temperature, whose convergence guaran- 
tees the existence of the bound state. However, definite 
conclusion on whether Xc melts cannot be drawn from 
the current consideration. 

Solving Eq. ([9]) with respect to m and F for the mini- 
mum of OPE side moment ratio at each temperature, we 
obtain the relation between mass shift and width broad- 
ening. The results are displayed in Fig. 21 As in the 
S'-wave cases, a linear relation is also seen in the P-wave 
quarkonia. In the small mass shift and large broadening 
case at r = l.OSTc, the relation shows some deviations 
from the linear relation. This comes from the fact that 
the stability is achieved at large n which gives strong 
dependence of the moment ratio on n. As mentioned be- 
fore, convergence of OPE is not so good. Strictly speak- 
ing, physical parameters extracted from the sum rules 
should not strongly depend on the external parameters, 
^ and n. In this case, however, F at Sm = depends 
on n while Sm does not. Therefore the results of the 
width might be less reliable. In reality there should be 
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FIG. 2: (Color online). OPE terms as a function of T/T^. At 
each temperatures, the terms are calculated at n which gives 
the minimum of the moment ratio as shown in Fig. [T] 
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FIG. 3: (Color online). OPE terms as a function of n for 
vacuum and T = l.OSTc. 



mass shifts in this temperature region caused by change 
of the chromo-electric field, i.e., second-order Stark ef- 
fect [13, H^. Hence, a small mass shift with large width 
broadening would not be a realistic combination. 

To see the temperature dependence, we depict the tem- 
perature dependence of the maximum mass shift and 
width broadening in Fig. [3 One can see that both Xco 
and Xci exhibit a clear critical behavior with respect to 
the abrupt change of the gluon condensates. This qual- 
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FIG. 4: (Color online). Relation between the mass shift and 
width broadening at temperatures above Tc- The upper and 
lower panels show XcO and Xci cases, respectively. 



itative feature was also seen in the S'-wave [i2[ cases, 
but in the present case, the spectral changes are roughly 
a factor of two larger . Especially the maximum mass 
shifts reach -150 MeV at T = Tc and more than 400 
MeV at r = l.OSTc- Lattice calculation suggests [8] that 
the XcO a-iid Xci will dissolve before T = l.lTc- However, 
as is evident also in the lattice data, our result suggest 
that just before it dissolves, the spectral density will be 
greatly modified and the mass greatly reduced slightly 
above which are more promising signals for direct ob- 
servation in future experiments than the correspondingly 
smaller mass shift expected for J/i/). 

We also calculate the mass shift and width broadening 
in nuclear matter. This can be done by replacing the 
change of the gluon condensates with those for the nu- 
clear matter [ij, . The result of mass shift and width 
broadening is shown in Fig. [S] The maximum mass shift 
is around 17 MeV, which is more than two times larger 
than that obtained previously for the S'-wave charmonia. 
The difference between XcO ^-nd Xci is around 1 MeV. 



IV. SUMMARY AND DISCUSSIONS 

In this paper, we extend the previous QCD sum rules 
analyses for the S'-wave charmonia to P-waves. The re- 
sults are qualitatively similar to those of S-waves but 
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FIG. 5: (Color online). Upper: mass shift as a function of 
T/Tc. Lower: Width broadening as a function of T/Tc. Cir- 
cles and triangles denote the XcO and Xci, respectively. 
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FIG. 6: Relation between mass shift and width broadening in 
the nuclear matter. Solid and dashed lines stand for XcO and 
Xci, respectively. 



show larger spectral changes. Due to the absence of con- 
vergence of the OPE, our approach is limited to tem- 
peratures above but close to Tc. Note that our results 
do not necessarily imply the absence of the bound states; 
spectral function obtained in lattice calculations does not 
have a good resolution of the width. Although potential 
models are extensively studied with potentials extracted 



from lattice calculations [2a, |27|, [23, [29[ , it is not straight- 
forward to implement the critical change into potentials 
at fixed temperatures, and also ambiguities exist in the 
choice of the potential. 

One may wonder if the continuum part in the hadronic 
spectral function that we have neglect in this paper can 
affect the result. As described before, contribution from 
this part can be suppressed by making use of the moment 
sum rule prescription. However, this suppression lead to 
too large n of the moment at T > l.OSTc, which results 
in the non-monotonic behavior of the width in Fig. [5j We 
will be able to improve it by taking the continuum into 
account. In this case, however, the continuum threshold 
So in the phenomenological side appears as a free pa- 
rameter. This can be fixed so that mass agrees with the 
experimental value in the vacuum case, while we have no 
criteria in the case of medium. This means the change of 
the OPE in medium might be attributed to the change 
of the threshold. Nevertheless, present calculation has 
general consequences because it corresponds to results 
at So ^ oo limit. Therefore, Fig. [T] shows the maxi- 
mum mass at each temperatures, as it is well known in 
the charmonium sum rule (20| that the extracted mass 
increases with increasing sq- Consequently, if one ad- 
justs So so as to reproduce the experimental value at 
zero temperature and then keep it fixed as temperature 
changes, the analysis will result in similar amount of spec- 
tral changes. If one allows Sq to change with tempera- 
ture, the mass shifts will be larger than those obtained 
at So — > oo. Even in this case, there are clear evidence of 
the change of the masses at T > 1.02Tc at which masses 
are smaller than experimental values. Further discussion 
on the effect of the continuum will be given in a future 
publication [30| . Future works also include directly com- 
paring the QCD sum rule result with the T-dependence 
of the Euclidean correlation functions observed in lattice 
QCD calculations. 

In this paper, we utilized the gluon condensates ex- 
tracted from lattice gauge theory of the pure gluonic 
system. Since there are many differences between the 
quenched and full QCD, one may wonder if the present 
results are relevant for realistic situations. First, the ef- 
fect of dynamical quarks alters thermodynamic proper- 
ties. Namely, transition temperature reduces to 170-200 
MeV [m while it is around 260 MeV for pure SU(3) H. 
On the other hand, what is important for the present 
analysis is the the quantitative change in magnitude of 
the condensate near Tc- Figure 21 of Ref. [T§| shows the 
behavior of the scalar gluon condensates of full QCD case 
after the fermionic contribution is subtracted out. One 
sees that the amount of decrease from the vacuum value 
is almost the same as in the pure gauge case while the 
change around Tc becomes rather moderate. Consider- 
ing the fact that the present lattice results have been 
performed on quite coarse lattice for which the hadron 
spectrum is still distorted and thus the thermodynamic 
quantities in the low temperature region are not prop- 
erly calculated the changes of gluon condensate could 
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be larger below Tc [s^. This means that the spectral 
changes set in at much lower temperature below due 
to the large number of degree of freedom in the hadronic 
phase and the amount of change at Tc is almost the same. 
Although not available at present, the twist-2 gluon part 
will also change a little as it is proportional to the entropy 
density of the gluons, which will not be effected greatly 
by the dynamical quarks. Hence, our present result could 
change a little but the main features will remain if dy- 
namical quarks are introduced. 

The results contain rich physical consequences relevant 
to recent and future experiments. First of all, a large 
fraction of J /ip comes from Xc- The in-medium modifi- 
cation of Xc will affect the final yields of J/i/"- If both 
J/V' and Xc can exist slightly above Tc, Xc with modified 
mass and width will decay into J /ip- This will give fur- 
ther enhancement of J /ijj production. Moreover, if the 
charmonia will be produced statistically [34] at the same 
hadronization temperature, the larger mass shift of Xc 
will lead to a larger enhancement of the Xc compared to 
J I tp and then will result in the change of the production 
ratio Nj/^/N^^ in nucleus-nucleus collisions. Alterna- 
tively, since the width of Xc can be large enough to decay 
inside medium, it might provide a chance to directly mea- 
sure the large mass shift. This is experimentally a chal- 



lenging subject, which requires the reconstruction of the 
J/^ + 7 invariant mass distribution with extremely fine 
resolution, but if successful will provide a direct evidence 
of the critical phenomenon in QCD. Our result for the 
nuclear matter can be a testing ground of this idea. The 
maximally expected mass shift around 15 MeV may not 
be large, but enough to be observable in the anti-proton 
project at FAIR [35]. This will be another challenge for 
future experiments. 

To conclude, we have studied the in-medium spectral 
change of Xco and y^i with the QCD sum rule approach 
developed in Ref. [ll| . Our results show that the change 
in these states are more prominent than previously stud- 
ied S'-wave states, and this result will provide another 
means to probe critical phenomena in QCD. 
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